Electronic and optical properties of carbon nanodisks and nanocones 
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A theoretical study of the electronic properties of nanodisks and nanocones is presented within the 
framework of a tight-binding scheme. The electronic densities of states and absorption coefficients 
are calculated for such structures with different sizes and topologies. A discrete position approxi- 
mation is used to describe the electronic states taking into account the effect of the overlap integral 
to first order. For small finite systems, both total and local densities of states depend sensitively on 
the number of atoms and characteristic geometry of the structures. Results for the local densities 
of charge reveal a finite charge distribution around some atoms at the apices and borders of the 
cone structures. For structures with more than 5000 atoms, the contribution to the total density of 
states near the Fermi level essentially comes from states localized at the edges. For other energies 
the average density of states exhibits similar features to the case of a graphene lattice. Results for 
the absorption spectra of nanocones show a peculiar dependence on the photon polarization in the 
infrared range for all investigated structures. 



I. INTRODUCTION 

Carbon nanocones (CNCs) have been observed by mi- 
croscopy techniques [TJ [2], showing an ordered atomic 
structure. Large progress has been made on synthe- 
sis, characterization and manipulation of CNCs and car- 
bon nanodisks (CNDs) [3HH]. From the technological 
point of view, applications such as microscopy probes 
and electron-emitter devices may be envisaged by con- 
sidering the electric current established through the cone 
apex when electric fields are applied [9 . 

There are different theoretical schemes to describe 
the electronic properties of cone-like structures. Mod- 
els based on the Dirac equation [TU1 [UJ give a conve- 
nient insight of properties in the long wavelength limit. 
However, for finite-size graphenes, the longest stationary 
wavelength occurs in the border and a correct descrip- 
tion of the states near the Fermi level is given in terms 
of edge states [12 |T3]. Ab initio models [HHToj are able 
to predict detailed features, but they are restricted to 
structures composed of a few hundred atoms due to their 
considerable computational costs. Calculations based on 
a single tt orbital are able to describe the relevant elec- 
tronic properties [17l HE]- I n that spirit, we calculate 
the electronic structure and optical spectra of CNDs and 
CNCs within a tight binding approach. CNC structured 
systems generated by pentagonal and heptagonal defects 
were previously studied using a Green function recursive 
method [THUTS]. It was shown that, for cones generated 
with an odd disclination number n w , it is not possible 
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to define A and B graphene sublattices. In this case, 
therefore, the electron-hole symmetry is broken. 

The total number Nq of carbon atoms in a cone struc- 
ture may be estimated by dividing the cone surface area 
by half of the hexagonal cell's surface, 

N c = [4tt/(3a/3)](1 - n w /6){r D /a CC ) 2 , (1) 

where the disclination number n w corresponds to the in- 
teger number of 7r/3 wedge-sections eliminated from the 
disk structure, and is the cone generatrix [see Figure 
[TJ. The nanocone disclination angle is given by n w ir/3. 
For example, for n w = 1 and ro = 1 /im, the CNC has 
ss 10 s atoms. By extracting an integer number n w of 7r/3 
sections from a carbon-disk [cf. Figure [l], it is possible to 
construct up to 5 different closed cones. For n w = 1, the 
cone angle is 29i = 112.9°, corresponding to the flattest 
possible cone. In this case, h/rc — 0.66 and h/rry = 0.55. 

In this work, finite-size systems (from two hundred up 
to five thousand atoms) are studied by performing di- 
rect diagonalizations of the stationary wave equation in 
the framework of a first-neighbor tight-binding approach. 
Each carbon atom has three nearest neighbours, except 
the border atoms for which dangling bonds are present. 
The overlap integral s is considered different from zero. 
As we will show later, this has important effects on the 
cone energy spectrum. 

Although some of the graphene electronic properties 
are present in the CNCs, deviations are always mani- 
fested as a consequence of the different atomic arrange- 
ments, the finite-size of the nanocones, and also the pos- 
sible point symmetry of the distinct cones. In the ab- 
sence of external fields, the calculated density of states 
(DOS) shows a peak at the Fermi energy and the local 
density of states (LDOS) shows that electron states are 
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Figure 1: (Color online) Pictorial view of (a) a carbon disk 
composed of six wedge-section of 7r/3 central angle, then (b) 
the remotion of — {1, . . . , 5} wedge-sections from the 

disk, and (c) by folding it is constructed a cone. Geometrical 
elements: generatrix m, height h c , base radius r c , and apex 
opening angle 29, where sin 6 = 1 — n w /6. 



localized at the cone base. On the other hand, the sym- 
metries observed in the LDOS at different energies allow 
a systematic description of the electronic structure and 
selection rules of optical transitions driven by polarized 
radiation. Unlike the nanodisk, the presence of topolog- 
ical disorders in nanocones involve a deviation from the 
electrical neutrality at the apex and at the zigzag edges. 



II. THEORETICAL FRAMEWORK 

In what follows, we present results for n w = {0, 1,2}, 
corresponding to CND and CNCs whose disclination an- 
gles are 60° and 120°. For those systems, the sp 2 hy- 
bridization may be neglected. The electronic wave func- 
tion may be written as 



N c 



(2) 



where the denote the atomic orbitals 2p at site Ri. 
Note that the overlapping between neighbouring orbitals 
prevents the set consisting of the ket to be an or- 
thogonal basis. Only in the ideal case of zero overlap 
s = 0, the coefficients C° in = Ef=iC"k°> might 
be considered equal to the discrete amplitude probability 
(7r° | <!') to find an electron at the j-th atom (described by 
the one electron state |^)). We use the s ^ basis, 
to construct the eigenvalue equation and the \tt^} base 
to calculate the properties related to discrete positions. 
Of course, to relate both bases it is required to know the 
(7r?|7Tj) projection. 

We define a Nc X Nq matrix A^ 1 ' relating the nearest 
neighbouring atomic sites 



Similarly, 



i 1 1 J 1, are n.n. 

0, otherwise 



A (o) fl, (i = j) 
10, otherwise 



(3) 



(4) 



The S overlap matrix elements are then given by 

S ij = {n i \n j ) = A^+sA^. (5) 

The hopping matrix elements of the tight-binding Hamil- 
tonian H = + V are 



tA 



(i) 



(6) 



where t is the hopping energy parameter. Assuming the 
eigenvalue equation H( at ') — e2 P \^j), the atomic ma- 
trix elements are 



(n\HM\* j )=e 2p 



A^+sAW 



and 



(iri\H\9) =e(7Ti|*) , l<i<N c 



(7) 



(8) 



The resulting equation system may be written as a 
generalized eigenvalue problem HC = eSC, where the 
column vector C contains the coefficient Cj , 



e 2p (A^ + sA^) + tA^ 



C = e 



c. 



(9) 

The general solution may be expressed in terms of the 
auxiliar variables (7(0) and e(0), which satisfy 

iA«(7(0) = £ (0)C(0). (10) 

As (7(0) also satisfies Eq. Q, we obtain 

e = [ £2p + (1 + fleap/t)e(0)]/[l + se(0)/t] . (11) 

The orthogonality condition for the electronic states 

(V k \$ l ) = C ki S C l = 5 kJ (12) 

implies that 



C = 



g(Q) 

C(0) f 5(7(0) 



(13) 



For the calculation of the DOS we use a Lorentzian dis- 
tribution 



N c N c 

DOS(e) = 2^S(s- J) « 2 £ p2 

3=1 3=1 V ' 



T/tt 



(14) 



The LDOS is calculated in terms of the discrete am- 
plitude probability, (7r?|* 3 '), 

N c 

WOS(R i ,s) = 2j^\{^ j )\ 2 8(e-^) , (15) 

3=1 



where 



A (o) + £ A (i) 
2 



C 3 



(16) 
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as it is shown in Appendix A. 

The local density of charge (LDOC) related to the w 
electrons is calculated by assuming that the other 5 elec- 
trons and the 6 protons of the carbon atom act as a net 
charge +e. Assuming zero temperature and the indepen- 
dent electron approximation, only the states 1 < j < np 
will be occupied, where 



Up 



A c /2 , A c is even 
(N c + l)/2 , A c is odd. 



(17) 



Taking into account that the states below np con- 
tribute with — 2e and the fact that the np state con- 
tribution depends on the parity of the number of atoms 
in the system, the LDOC is written as 



LDOC(i? 4 ) = 



Tip — 1 



1-2 £ |(7r0|^)| 2 - 7 |(^|^)|2 

3=1 

(18) 

with 7=0 and 1, for Ac even and odd, respectively. 

Optical absorption coefficients a e (ui) are calculated by 
considering perpendicular (el = e x ,e y ) and parallel (eji = 
e z ) polarizations, in relation to the cone axis, 

ag(cj)oc -^2\{V l \e-p\W)\ 2 d^ 1 -e l -hu] , (19) 



with e l,J corresponding to the energies of occupied and 
unoccupied states, respectively. 

The oscillator strength may be written in terms of the 
spatial operators (x, y and z) |20j . i.e., 



(¥\p z \W) = ^{e> - e*)(W\z\V) , (20) 
in 

where (^> l \z \W) is calculated to first order in s, using 
(|A9| of Appendix A, 



(mm 3 ) 



c 



*t 



: + - ( A (1) z + zA (1) 



C 3 



(21) 



III. RESULTS AND DISCUSSION 

A. Electronic Density of States 

In what follows, we present numerical results for sys- 
tems composed of up to 5000 atoms. In the limiting 
case of Ac —> oo, the energy spectrum is in the range 
from e min = -3|i|/(l + 3s) to e max = +3|i|/(l - 3s), 
the van Hove singularities occur at = — \t\/(l + s), 
£^ H = +|t|/(l — s), and the Fermi energy is at ep = 0. 
A r = 1 1 1 /L00 broadening and an overlap s = 0.13 are 
assumed. Different plots in Figure [2] show the density 
of states averaged over the Ac atoms and the LDOS for 
a CND [Figures [2Ta) and[2jd)], a single-pentagon CNC 
[Figures [2jb) andMe)], and for a 2-pentagon CNC [Fig- 
ures ^c) and|2|;f)]7for A c = 258,245, and 246, respec- 
tively. All results are shown in an energy range around 
£2p = 0. 
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Figure 2: (color online) DOS and LDOS for a N c = 258 
nanodisk [(a) and (d)], a Nc = 245 1-pentagon nanocone [(b) 
and (e)], and a Nc = 246 2-pentagon nanocone [(c) and (f)]. 
LDOS curves (black-1, red-2, blue-3) match, for each system, 
the corresponding atoms shown in Figure [3] Vertical lines in 
each panel indicate the position of the Fermi energy. 




Figure 3: Pictorial view of (a) a nanodisk center, (b) a 
1-pentagon nanocone apex, and (c) a 2-pentagon nanocone 
apex. Atoms with different colors-numbers indicate different 
point symmetries for each system. 




Figure 4: (color online) DOS and LDOS for a N c = 5016 
nanodisk [(a) and (d)], a Nc = 5005 1-pentagon nanocone 
[(b) and (e)], and a Nc = 5002 2-pentagon nanocone [(c) and 
(f)]. LDOS curves (black-1, red-2, blue-3) match, for each 
system, the corresponding atoms shown in Figure [3] Vertical 
lines in each panel indicate the position of the Fermi energy. 



4 




Figure 5: (Color online) LDOS in arbitrary units for a 5016-atom nanodisk [panels (a) to (e)], a 5005-atom nanocone with one 
pentagon at the apex [panels (f) to (j)], and a 5002-atom nanocone with 2 pentagons at apex [panels (k) to (o)]. The considered 
energies are (a,f,k) e m i n , (b,g,l) e" H , (c,h,m) ep, (d,i,n) e£ H , and (e,j,o) e max . The LDOS is measured with respect to the mean 
LDOS which is equal to the DOS at the considered energy. 



As expected for small finite systems, the DOS, LDOS, 
and the position of the Fermi energy depend on the num- 
ber of atoms considered in the numerical calculation and 
on their characteristic geometries [2"TM23| . A remarkable 
difference between CND and CNCs structures is the exis- 
tence of a finite DOS above the Fermi level for nanocones. 
This clear metallic character of the DOS for nanocones 
is more robust for the 2-pentagon CNC PJU [21] • One 
may remark that this feature comes out from a symme- 
try break induced by the presence of topological defects 
in the CNC lattices, which generates new states above 
the Fermi energy, not present in the CND structure. The 
contributions to the DOS coming from the apex atoms 
states are apparent in the LDOS of Figures [2j[e) andFJJf). 
Also notice that for the 2-pentagon case, in which there 
is a large topological disorder, the LDOS spectra exhibit 
significant differences depending on the point symmetry 
of the considered atom (cf . Figure [3| . 

For increasing number of atoms, the total DOS for 
the different nanostructures are very similar to the corre- 
sponding DOS of a graphene layer, except for the edges 
states which show up as a peak at the Fermi energy, as 
shown in Figures (4^a), |4^b), and|4^c). ft is interesting 
to note that the apex atomic states do not contribute to 
the total DOS near the Fermi energy but mainly near the 
graphene- like van Hove peaks. Notice that in the case of 
two-pentagon nanocones the LDOS at the tip exhibits a 
robust metallic character. 

To analyse the finite-size effects and the role played by 
the different symmetries of the cone-tip sites, we depict 
LDOS contour plots for the three studied structures by 
considering some characteristic energies: the minimum 
energy, the resonant peak below the Fermi energy, the 



Fermi energy, the resonant peak above the Fermi energy, 
and the maximum energy. Figure [5] illustrates the ex- 
ample of a CND with 50f6 atoms (top row), a single- 
pentagon CNC with 5005 atoms (middle row), and a 
2-pentagon CNC with 5002 atoms (bottom row). The 
electronic states corresponding to energies at the band 
extrema have the largest wavelength compared to the 
characteristic size of the system. In this way, the de- 
tails of the lattice become less important and the states 
exhibit azimuthal symmetry. An interesting feature for 
the nanocones is that at these energies the apex corre- 
sponds to a node for the maximum energy and an antin- 
ode for the minimum energy, respectively. On the other 
hand, the states at the Fermi energy are localized at the 
cone border, mainly at the zigzag edges as it is clearly 
shown in Figures [5^c), [5jh), and[5]jm). For the states 
whose energy is near to the van Hove peaks, the LDOS 
reflect the symmetries of each system, i.e., for CND the 
27r/6-rotation symmetry and 12 specular planes [cf. Fig- 
ures U» andgd)], for a single-pentagon CNC there is a 
27r/5-rotation symmetry and 5 specular planes [cf. Fig- 
ures [5jg) and[5ji)], and for a two-pentagon CNC, there 
is a 7r/2 rotation symmetry and 2 specular planes [cf. 
Figures [BJl) and[5ji)]. 



B. Electric Charge Distribution 

The electric charge per site, in terms of the fundamen- 
tal charge e, was obtained using Eq. (18 1. Results for 



the electric charge distribution for CNDs indicate that 
all the atomic sites preserve the charge neutrality, i.e., 
LDOC=0. For the CNCs, however, the atoms at the 
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sites 


1 


2 


3 


max 


1-pentagon 


-0.071e 


+0.014e 


-0.059e 


+0.042e 


2-pentagon 


-0.055e 


-0.067c 


-0.066e 


+0.076e 



Table I: LDOC (fundamental charge units) at some relevant 
atoms in the cone apices shown in Figure [gb) and (c). The 
maximum (max) value occurs at the zigzag edge of each sys- 
tem. 




Figure 6: (Color online) Electric charge distribution in neutral 
CNCs, for a single-pentagon structure with 245 atoms (a) and 
two-pentagon system with 246 atoms (b). 



apex acquire negative charge and the atoms around the 
cone base exhibit positive charges at the zigzag edges. 
As Nc increases, the LDOCs at the apices, for the two 
studied CNC structures, tend to the asymptotic values 
shown in table [TJ which are in good agreement with the 
values reported by Green method calculations [T51 HH] . 

Figure [6] depicts the LDOC for the two types of CNC 
structures, showing that the nonequilibrium of the charge 
distribution is restricted to the apex and edge regions: 
electric neutrality is found at all the other surface sites. 
The values found for the LDOC at the apex regions are 
found to be independent of the size of the cones whereas 
this is not true for the edge states. When the number 
of atoms of the CNC structure is even, the edge-states 
LDOC exhibits the same symmetry of the cone. For odd 
Nc, the Fermi level is occupied by a single electron and 
then the LDOC at the edge states reflect the broken sym- 
metry. 



C. Absorption Spectra 

We have also calculated the absorption coefficient for 
the CND and CNC structures, for different photon po- 
larizations. Figure [7] shows the results for the absorption 
coefficients a x and a y , for polarization perpendicular to 
the cone axis, and a z for parallel polarization. Calcu- 
lated results are shown for a nanodisk composed of 5016 
atoms, a single-pentagon nanocone with 5005 atoms, and 
a 2-pentagon nanocone with 5002 atoms. For the case of 
large CNDs, the spectra present the general features ob- 
served for the absorption of a graphene monolayer. In 
the infrared region the absorption coefficient of a gra- 
phene monolayer is expected to be strictly constant [25], 
whereas for higher energies the spectrum shows a strong 
interband absorption peak coming from transitions near 
the M point of the Brillouin zone of graphene [2B]. The 




0.1 

o.o 

0.3 
?0.2 



j, «*"«y 


■ 1. 


i 

<b) 


1.1,1.1,1. 






CO 





















1.0 1.6 2.0 

photon energy 



Figure 7: (Color online) Absorption coefficient for x (black 
curves), y (red curves) and z (blue lines) polarizations for (a) 
a nanodisk with 5016 atoms, (b) a single-pentagon nanocone 
composed of 5005 atoms, and (c) a two-pentagons nanocone 
with 5002 atoms. The photon energies are given in units of 

hoj/t. 



main difference for a finite CND is a departure from a 
completely frequency-independent behavior for low ener- 
gies, where the absorption coefficient shows oscillations 
as a function of the photon energy instead of a constant 
value. This is a consequence of the border states that 
are manifested as a peak in the total DOS at the Fermi 
energy [371 HI] ■ F° r CNCs, the general behaviour is the 
same as for nanodisks, except for the dependence of the 
absorption on the photon polarization, in particular for 
low energies. Furthermore, the main absorption peaks 
for different polarizations occur when the photon energy 
is equal to the energy between the two DOS van Hove- 
like peaks (cf. Figure[4|. Notice that the overlap integral 
s ^ leads to an energy shift of the main resonant ab- 
sorption peak given by 5 « 2s 2 |t|/(l — s 2 ) w 100 meV. 
This is a significant value for actual experimental mea- 
surements. 

Concerning the different polarization directions, one 
should notice that, as occurs in Cq v symmetric systems, 
a z = and a x — a y for the nanodisk. On the other 
hand, the absorption coefficients for the different cones 
studied (single and two pentagons) are finite for parallel 
polarization, and it depends on the structure details: as 
a z increases for a two-pentagon CNC structure, ol x v de- 
creases. Due to the lack of n/2- rotation symmetry, one 
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Figure 8: (Color online) Absorption coefficient for x (black 
curves), y (red curves) and z (blue lines) polarizations for (a) 
a nanodisk with 258 atoms, (b) a single-pentagon nanocone 
composed of 245 atoms, and (c) a two-pentagon nanocone 
with 246 atoms. The photon energies are given in units of 
hijj/t. Curves in each panel are vertically shifted, for better 
visualization of different polarization results. 



should expect, in principle, different results for x- and 
y-polarizations for any nanocone. However, such differ- 
ence is observable just for the absorption coefficient of 
the two-pentagon CNC system, mainly in the range of 
low photon energies. The fact that a x = a y , for the case 
of 1-pentagon CNC structure, may be explained using 
similar symmetry arguments applied to Cq v symmetry 
dots [28], extended to the C$ v symmetric cones. In the 
case of a 2-pentagon CNC, the apex exibits a Ci v sym- 
metry, preventing the cone to be a C^ v symmetric sys- 
tem. As the apex plays a minor role, a x and a y will be 
slightly different. A large difference between the a z and 
the a x ,y CNC absorption spectra occurs in the limit of 
low radiation energy. The a z coefficient goes to zero as 
hoj — s- whereas a x<y shows oscillatory features. The be- 
haviour of the absorption for parallel polarization is due 
to the localization of the electronic states at the atomic 
sites around the cone border. As the spatial distribution 
of those states are restricted to a narrow extension along 
the z coordinate, the z degree of freedom is frozen for low 
excitation energies. 

The dependence of the absorption spectra on the ge- 
ometrical details of the different structures is more no- 
ticeable for finite-size nanostructures. This can be seen 



in Figure [8] which depicts the absorption coefficients for 
the CND composed of 258 atoms, the single-pentagon 
CNC with 245 atoms, and the 2-pentagon CNC with 
246 atoms. The degeneracy of the x- and y- polariza- 
tion spectra is apparent for the smaller one-pentagon na- 
nocone, as expected due to symmetry issues. On the 
other hand, the symmetry reduction for the 2-pentagon 
structure leads to a rich absorption spectra, exhibiting 
peaks at different energies and with comparable weights 
for distinct polarizations. In that sense, one may pro- 
pose absorption experiments as an alternative route to 
distinguish between different nanocone geometries. 



IV. CONCLUSIONS 

Here we have presented a theoretical study on the 
electronic properties of nanodisks and nanocones in the 
framework of a tight-binding approach. We have pro- 
posed a discrete position approximation to describe the 
electronic states which takes into account the effect of 
the overlap integral to first order. While the \tt) base 
keeps the phenomenology of the overlap between neigh- 
boring atomic orbitals, the \tt°) base allows the construc- 
tion of diagonal matrices of position-dependent opera- 
tors. A transformation rule was set up to take advantage 
of these two bases scenarios. 

We have investigated the effects on the DOS and 
LDOS, of the size and topology of CND and CNC stru- 
tures. We have found that both total and local density 
of states sensitively depend on the number of atoms and 
characteristic geometry of the structures. One impor- 
tant aspect is the fact that cone and disk borders play 
a relevant role on the LDOS at the Fermi energy. For 
small finite systems the presence of states localized in the 
cone apices determines the form of the DOS close to the 
Fermi energy. The observed features indicate that small 
nanocones could present good field-emission properties. 
This is corroborated by the calculation of the LDOC, 
that indicates the existence of finite charges at the apex 
region of the nanocones. For large systems, the contri- 
bution to the DOS near the Fermi level is mainly due to 
states localized in the edges of the structures wheres for 
other energies, the DOS exhibits similar features to the 
case of a graphene lattice. 

The absorption coefficient for different CNC structures 
are calculated and we have found a peculiar dependence 
on the photon polarization in the infrared range for the 
investigated systems. The symmetry reduction of the 2- 
pentagon nanocones causes the formation of highly struc- 
tured absorption spectra, with comparable weights for 
distinct polarizations. The breaking of the degeneracy for 
different polarizations is found to be more pronounced for 
small nanocones. Absorption experiments may be used 
as natural measurements to distinguish between different 
nanocone geometries. 
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Appendix A 

A discrete position scheme in terms of the |-7T°) states 
was used to represent functions of the position given in 
terms of the atomic base, since they satisfy the same 
properties of the position states, i.e., orthogonality 



,(0) 



and completeness 



N c 

fe=l 



(Al) 



(A2) 



in a Ac-dimensional subspace. The identity operator 
may also be constructed using the s ^ base as 



i = ^K)(5- 1 ) fe , 



kl TO | , 



(A3) 



k,l 



with the S' 1 « A^ ) - sA^ + 0{s 2 ) matrix being dif- 
ferent from the Nc x Nc identity matrix A^ ^ . 

We take \n°) as the discrete position state and assume 
that the matrix elements f3 of position-dependent func- 



/^ = K°l/V°)=/(4)Ag ) . (A4) 

Differently from the f R matrices, / matrices in the s ^ 
representation 



fij = (lil/Ki) 



(A5) 



are not diagonal. However, by performing the similarity 
transformation 

fol/Vi) = <^l/h°> (T?ki> > (A6) 

k,l 

we may obtain the unknown / matrix in terms of the 
known f R matrix, provided the transformation rule be- 
tween the 7T° and ir bases is known. By assuming 
(7if |7Tj) = aA\j + /3A-j\ the s ^ representation may 
be found. The coefficients a and j3 are obtained by using 
the identity |A2) into Eq. @, 



N c 



k=l 

and, to first order in s, (a = 1 and /3 = s/2) we have 

(7r°K-) = AW + ( s /2)A«. (A8) 



By replacing (|A8|) in (|A6|), one obtains 

Hi 



f (R) + £ / /( H) A (1) + A (l) + 0(s 2 } 
2 V / J 



(A9) 

as the matrix elements of a position-dependent function 
in the 7r-base. 
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